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(N ■ Abstract 
> 
On 
O . 

, We investigate the dynamics of two point-like particles through the third 

post-Newtonian (3PN) approximation of general relativity. The infinite self- 

o 



field of each point-mass is regularized by means of Hadamard's concept of 
"partie finie" . Distributional forms associated with the regularization are used 
systematically in the computation. We determine the stress-energy tensor of 



point-like particles compatible with the previous regularization. The Einstein 
field equations in harmonic coordinates are iterated to the 3PN order. The 
3PN equations of motion are Lorentz-invariant and admit a conserved energy 
(neglecting the 2.5PN radiation reaction). They depend on an undetermined 
coefficient, in agreement with an earlier result of Jaranowski and Schafer. This 
suggests an incompleteness of the formalism (in this stage of development) 
at the 3PN order. In this paper we present the equations of motion in the 
center-of-mass frame and in the case of circular orbits. 



1 



Typeset using REVT^X 



A cardinal problem in Gravitational Physics is that of the dynamics of binary systems 
of point particles. In general relativity, this problem is tackled by means of the post- 
Newtonian approximation, or formal expansion when the speed of light c goes to infinity. 
By definition, the nPN approximation refers to the terms in the equations of motion that 
are smaller than the Newtonian force by a factor of order l/c 2n . For the motion of two non- 
spinning point particles, the 1PN approximation was obtained first by Lorentz and Droste 
Hfl. Subsequently, Einstein, Infeld and Hoffmann f2| re-derived the 1PN order using their 
famous "surface-integral" method. In the eighties, Damour and Deruelle |§, starting from 
a "post-Minkowskian" iteration scheme developed by Bel et al , were able to compute the 
equations of motion up to the 2.5PN order, at which the gravitational-radiation reaction 
effects first take place. The motivation was to firmly establishing the rate at which the 
orbit of the binary pulsar PSR 1913+16 decays because of gravitational-radiation emission. 
The 2.5PN approximation was then obtained by Schafer || using an "ADM Hamiltonian" 
approach initiated by Ohta et al ||. Furthermore, Kopeikin et al |7| derived the same 
result within their "extended-body" method (without any need of a regularization) . More 
recently, the 2.5PN equations of motion as well as 2.5PN gravitational field were derived 
by Blanchet, Faye and Ponsot [|5] applying a direct "post-Newtonian" iteration of the field 
equations. Finally, Jaranowski and Schafer || investigated within the Hamiltonian approach 
the 3PN order and found some ambiguities linked to the regularization of the self-field of 



point masses. As for them, the 3.5PN terms in the equations of motion are well-known pO 
they are associated with higher-order radiation reaction effects. 

The motivation for working out the 3PN equations of motion is not the timing of the 
binary pulsar anymore, but the detection of gravitational radiation by future experiments 
such as LIGO and VIRGO. Indeed, the 3PN equations are needed in particular to construct 
accurate 3.5PN templates for detecting and analyzing the waves generated by inspiralling 
compact binaries |IT[ . Currently, we know the complete templates up to the 2.5PN order 
T2] |, plus the contribution therein of non-linear effects to the 3.5PN order [ T3(] , plus all the 



terms in the vanishing mass-ratio limit to the 5.5PN order ||14|| . In this Letter, we outline our 



derivation of the 3PN equations of motion, based on the direct post-Newtonian approach 
of ||, and we present the result in the case, appropriate to inspiralling compact binaries, 
of circular orbits. We confirm by means of a well-defined regularization a la Hadamard the 
finding of Jaranowski and Schafer || that there remains at the 3PN order an undetermined 
coefficient appearing in front of a quartically non-linear term (proportional to G 4 ). 

Consider the class T of functions F(x) that are smooth on M 3 except at two isolated 
points yi and y2, around which they admit a power-like singular expansion of the form 

VnGN, F(x) = r i /a(ni) + o(r") when r x -»• , (1) 

a,Q<a<n 

where t\ = |x — yi|, ni = (x — yi)/rj, and where the powers a are supposed to be real, to 
range in discrete steps: a G (aj)igN, an d to be bounded from below: a < a - The coefficients 
if a of the various powers of r x in this expansion are smooth functions of the unit vector n^. 
We refer to the coefficients if a with a < as the singular coefficients of F around 1; their 
number is always finite. Moreover, we have the same type of expansion around the other 
point (when r 2 — > 0). The Hadamard "partie finie" [15] of F at the location of the singular 
point 1 is equal to the angular average of the zeroth-order coefficient in (|l]), i.e. 

M = J f£ /»(".) • (2) 

with dQ x = dQ(ni) the usual solid angle element. The partie finie is "non-distributive" in 
the sense that (FG)i ^ (F)i(G)i in general. Besides we define also the partie finie 
(Pf) of the divergent integral f o? 3 x F, assuming that F decreases sufficiently rapidly when 
|x| — > +oo so that the divergencies come only from the singular points 1 and 2. With full 
generality |15| , [T6H , 



Pf Sl)S2 d 3 xF = lim \ I rf 3 x F 

+ ^E o £^(9 1 +4 ' ln fe) (r ' F) ' + 1 " 2 }' (3) 

The first term is the finite integral over R 3 deprived from the two spherical balls Bi(s) and 
#2(5) with radius s and centred on the two singularities. The extra terms are such that 



they exactly cancel out the divergent part of the integral when s — > (the notation 1 <-> 2 
indicates the same extra terms but referring to the other singularity point). The logarithmic 
terms depend on two arbitrary positive constants s\ and S2 that come from the arbitrariness 
in the choice of unit length for measuring s; hence, the partie finie depends on both si and 
s 2 (as indicated by the notation Pf sl)S2 ). Applying (|3|) to the case of a gradient c^F, we find 



In words, the integral of a gradient is equal to the sum of the surface integrals surrounding 
the two singularities, in the limit where the surface areas shrink to zero and following the 
regularization (0). Thus, the integral of a gradient is not zero in general, which shows that 
the "ordinary" derivative <9,F is not adequate for applying to point-particles a formalism 
initially valid for continuous sources, since in the latter case the integral of a gradient 
does never contribute. To define a "better" notion of derivative, we must construct the 
distributional forms associated with the functions in the class T . 

For any F G F, we consider the "pseudo-function" PfF defined as the linear form on F 
such that VG G F, < PfF, G >= Pf f <i 3 x FG, the duality bracket denoting here the result 
of the action of PfF on the function G. The product of pseudo-functions is defined to be 
the ordinary pointwise product that we use in Physics, i.e. PfF.PfG = Pf (FG). With the 
help of the Riesz [[HJ delta-function <5 e (x — y x ) = £ ^~ £ ^ rf~ 3 , which belongs to the class F, 
we construct Jl8| the pseudo-function Pf<5i (in the limit e — > 0); by definition: VF G F, 
< PfSi,F >= (F) 1 . Clearly Pf^ generalizes the standard Dirac distribution 5i = <5(x — y x ) 
to the case of the Hadamard regularization of the functions in T . Furthermore, consistently 
with the product of pseudo-functions, we construct the object Pf(F^i) which is such that 
VG G J 7 , < Pf(FSi),G >= (FG)i. A trivial consequence of the non-distributivity of the 
Hadamard partie finie is that Pf (F5\) ^ (F)iPf5i in general cases. The derivative of the 
pseudo-function PfF is then obtained from the requirements that (i) the "rule of integration 
by parts" is satisfied, i.e. VF, G G F, < ^(PfF),G >= - < ^(PfG),F >, (ii) the 
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derivative reduces to the "ordinary" one in the case where all the singular coefficients of F 
vanish. These requirements imply in particular that < <9j(PhF), 1 >= 0, i.e. the integral of 
a gradient is zero. A derivative operator satisfying (i) and (ii) is given by ]Tj| 

9i(PfF) = Pf (V + 4vrnl k n f_ x + J] 1 /_ 2 _J ^ + 1 <- 2) (5) 
\ Li y 1 -] 1 J / 

fc>0 1 

(assuming for simplicity that the i/ 's have a G Z). This derivative reduces to the standard 
distributional derivative of Schwartz |L6 when applied on smooth functions with compact 



support. We refer to [18] for the construction of the most general derivative operator satisfy- 
ing (i), (ii) and, in addition, (iii) the rule of commutation of derivatives [not obeyed by (|5|)]. 
One can show however that it does not satisfy in general the Leibniz rule for the derivative 
of a product. The derivative (|5]) is sufficient in the derivation of the results below. See |18| 
for details about the Hadamard regularization and the associated pseudo-functions. 

In the post-Newtonian application we are led to consider the partie finie, in the sense 
of (|3p, of the Poisson integral of F, i.e. Pf J g? 3 xF(x)/|x — x'|; more specifically, we are 
interested in the regularized value, in the sense of (§), of the latter Poisson integral at the 
location of the singular point 1, i.e. when r[ = [3 



Pf 



S1,S2 



dh 



X — X 



Pf 
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F -4tt 



0. We obtain 18 



ln(' J 



r 

\Si 



4.1n(^)(^) 1 + 4.1n(^)(^) 2 + 



(6a) 
(6b) 

(with T12 = |yi — y2 1 ) • The first term in (|6aD represents simply what we get by replacing 
formally x' by y x inside the integrand of the Poisson integral. The second term is due to 
the presence of some logarithms lnr^ in the expansion of the integral. (An adaptation of 



the previous formalism, detailed in ||18|| , is needed to take these logarithms into account, as 
well as the presence of the integrable singularity x'.) As, at last, the lnr^ can be gauged 
away, we regard it as a constant, taking some finite value (even though r[ — > 0). We check, 
on the other hand, that the true constant Si cancels out between the two terms of (pa]), so 
that the result depends only on lnr^ and lns2- The complete dependence of the partie finie 
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on these constants is shown in (|6b"D, with the convention that the dots indicate the terms 
independent of the constants. 



The Einstein field equations relaxed by the condition of harmonic coordinates [i.e. 
d v bT = with bT = sf=gg» v - rf\ g = detg^; tf v =diag(-l, 1, 1, 1)] read as 



where T^ v is the matter stress-energy tensor and a complicated functional of h which 
is at least of order 0(h 2 ) (and where □ = rj^d^d^). We start by constructing a post- 
Newtonian solution of (0), initially valid in the case of a continuous ("fluid") matter tensor 
T^, and parametrized by some appropriate potentials. We define a "Newtonian" potential 
V = 4:nGa] where D^ 1 denotes the standard retarded integral and cr = (T 00 + T n )/c 2 ; 

we also introduce a 1PN "gravitomagnetic" potential Vi = 47rG<Tj] where <jj = T 0i /c; 

some 2PN potentials X, Ri and Wy, e.g. Wy = D^ 1 [— 4nG(o-ij — 5ij(Tkk) — diVdjV] where 
oij = T % i\ and finally some 3PN potentials T, and Zij. In particular, the potential 
generates the non-linear term n^ 1 [H/j J -9j J V A ], involving a cubic (G 3 ) contribution, which is 
part of the potential X (many other cubic terms are contained in T and Y\). With a specific 
choice of potentials we can arrange that all the quartic (G 4 ) terms in the metric appear in 
"all-integrated" form. Since V is dominantly Newtonian, it needs to be evaluated at the 
3PN order but, for instance, the term D^fW'ijdy V], inside the 2PN potential X, needs only 
a relative 1PN precision. The metric is expressed as a functional of all these potentials; and 
with our particular choice of potentials, it turns out not to be too complicated. 

An important point is now to determine the expression of the matter stress-energy tensor 
T^ v appropriate to the description of point-particles. We demand that the dynamics of 
point-masses follows from the variation, with respect to the metric, of the action 



where = (c,dyi/dt) is the coordinate velocity of particle 1. We can check that to the 
3PN order all the metric coefficients g^ v belong to T (treating lnr^ as a constant); so {g^ v )\ 



uhr = i^G(_^ r ^ + a^"[/i, dh, d 2 h] 



(7) 
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in (§) denotes the value of the metric at 1 in the sense of (0) [or, rather, in the sense of a 
Lorentz-covariant Hadamard regularization defined below]. The stationarity of the action 
with respect to a metric variation within the class T (i.e. Sg^ G T) yields the stress-energy 
tensor 

point - particle " v^i^MP \v=9J 

where the pseudo-function Pf(- 7 =5 1 ) is of the type Pf(F5i) defined before. [From the rule 
of multiplication of pseudo-functions we find that the matter source term in (0) involves 
the pseudo-function Pf (y/— go~i)-] To obtain the equations of motion of the particle 1 we 
integrate the matter equations of motion Vi/T^ nt _ particle = over a volume surrounding 1 
(exclusively), and use the properties of pseudo-functions. The equations turn out to have 
the same form as the geodesic equations, not with respect to some smooth background but 
with respect to the regularized metric generated by the two bodies. Namely, 

dt W-iap.Mv? J 2 V-ia^iv^ ' 

where all the quantities at 1 are evaluated using the regularization. Let us emphasize that 
the equations of motion ( [TDD are derived from the specific expression (^) of the stress-energy 
tensor; had we used another expression, e.g. by replacing Pf(- 7 =5 1 ) — > (- 7 =) 1 Pf5 1 inside 
(§) (which is forbidden by the non-distributivity of Hadamard's partie finie), we would have 
obtained some different-looking, and a priori uncorrect, equations. 

The regularization is defined stricto sensu within the spatial slice t =const, and 
therefore should prevent, at some stage, the equations of motion from being Lorentz invariant 
(recall that the harmonic-gauge condition preserves the Lorentz invariance). It is known 
that regularizing within the slice t =const yields the correct, Lorentz-invariant, equations 
of motion up to the 2PN level |J. We find that the breakdown of the Lorentz invariance 
due to the regularization (0) occurs precisely at the 3PN order. Therefore, starting at this 
order, we must in fact apply a Lorentz-covariant regularization in (PD-(|T0]) . Evidently, the 
good thing to do is to apply the Hadamard regularization in the frame at which the particle 



is instantaneously at rest. Let us consider the Lorentz boost x ,fJ- = A fJi u (V)x u , where V 
denotes the constant boost velocity. We replace all the quantities in the original frame 
by their equivalent expressions, developed to the 3PN order, in the new frame. Notably, 
the trajectories yi(i), y2(t) and velocities vi(t), v 2 (t) are replaced by certain functionals 
of x' and the new trajectories y^it'), y' 2 (t') and velocities v[(t'), v 2 (i') (where t' = x'°/c). 
We apply the Hadamard regularization within the slice t' =const, keeping V as a constant 
"spectator" vector. Finally, we re-express all the quantities back into the original frame at 
the point 1 (ri — > 0), under the condition that v[(t') = and (equivalently) V = vi(£). 
This ensures that the new frame is indeed the rest frame of the particule 1 at the instant t. 
Now the 3PN equations of motion are Lorentz invariant. 

All the potentials V, V i: Wij, ■ ■ ■ and their gradients are computed at the point 1, using 
the regularization of Poisson-type integrals defined by formulas like fl6a|). All the derivatives 
appearing inside the non-linear sources of the potentials are considered as distributional 
and evaluated following the prescription ([5]). We carefully take into account the fact that 
the distributional derivative does not obey the Leibniz rule (it does satisfy it only in an 
"integrated" sense, thanks to the rule of integration by parts). An important feature of 
the equations at the 3PN order is the occurence of some logarithms. From floTiD we know 
that they are necessarily of the type ln(ri2/r' 1 ) and In(ri2/s2) in the equations of motion of 
body 1; interestingly, the ln(r 12 /s2) appears only in a quartic-interaction term proportional 
to G^mim^. Thus, at this stage, the 3PN equations of 1 depend on the constants lnr^ and 
lns 2 (and idem for the equations of 2). Under the form we obtain them, the equations do 
not yet admit a conserved energy (of course we are speaking only about the conservative 
part of the acceleration, which excludes the radiation-reaction potential at 2.5PN order). 
However, we find that a conserved energy exists if and only if the logarithmic ratios ln(r 2 / s 2 ) 
and ln^/si) are adjusted in such a way that 



where A is a single numerical constant. If (and only if) the condition ( |lT|) is realized, the 




(11) 



S 



equations admit an energy and, in fact, a Lagrangian formulation; in this case, they depend 
on some arbitrary constant A. The dependence upon the masses in ([11]) is a priori allowed. 
Therefore, the formalism introduces at this point an undetermined constant A. [Using (|i~l~D, 
the equations of motion depend also on the constants lnr^ and lnr 2 , but it can be checked 
that the latter dependence is pure gauge.] 

Finally, having in view the application to inspiralling compact binaries, we present the 
3PN relative acceleration and center-of-mass energy in the case of circular orbits. The 
acceleration reads as 

^_. W _l F ,,„ + o(I), ( 12 ) 

where y i2 = yi — y 2 is the relative separation in harmonic coordinates, v 12 = dy^/dt 
the relative velocity, and F reac = _^ G3 ™ 3i/ Vl2 the standard radiation- reaction force at the 

r 12 

2.5PN order. The mass parameters are m = m\ + m 2 , = mim2 and v = — . The content of 



the 3PN approximation in ([12]) lies in the relation between the orbital frequency u and the 
coordinate distance r 12 = |yi 2 |. With 7 = -^3- denoting a small post-Newtonian parameter, 
we get 

Gm 



LU 2 



12 



|l + (-3 + z/) 7 + (6 + f z/ + z/ 2 ) 7 2 (13) 

22 In (^) +f A] z, + fz, 2 + z/) 7 3 + o( 7 4 )} . 



+ -10 + 



67759 _j_ 41 2 , 
840 ' 64 

The logarithm at 3PN depends on a constant r' Q defined as the "logarithmic" barycenter 
of the two constants r[ and r' 2 , namely lnr = — lnr^ + ^lnr' 2 . The constant r' can be 
eliminated by a change of coordinates. The center-of-mass energy E of the particles, such 
that ^ = as a consequence of the conservative equations of motion (neglecting F reac ), is 
obtained as 

E = ->^{l + (-1 + H 7 + (-1 + + X) 7 2 (14) 

+ (- w + [nr - w- 2 + f 1 - (t ) - f >\>> + 1- 2 + ii- 3 ) -r 3 + o(7*) } . 

At last, by substituting the expression of 7 in terms of the orbital frequency uo following 
from the inverse of ([13|), we find that the 3PN energy in invariant form is given by 

9 



E 



I /UC 2 x{l + (-§ 



3 _ ± u - 

4 12 



•)x+(-f + > 



J_ 

24 



2\ 2 
1/ J X 



(15) 



+ ( 



675 
64 



+ [ 



209323 
4032 



205^2 110 



96 9 



A] V 




■ 3 ) x 3 + 0(x 4 )}, 



with x = (^r^) 2 / 3 . In the form fll5j) the logarithm ln(r 12 /r ) cancels out. We can compare 
directly this result with the one obtained by Jaranowski and Schafer || (see the equation 



and osmetic = 24' wnere ^static and ^kinetic are the two "ambiguous" parameters found by 
Jaranowski and Schafer. Thus, our undetermined constant A defined by (|TTD is related to 
the ambiguous parameter cj sta ti C , while the other ambiguity cosmetic takes a unique value 



sense of Jaranowski and Schafer. Rather, the formalism is well-defined thanks in particular 
to the rules we employ for handling the pseudo-functions associated with the Hadamard 
regularization All the integrals encountered in the problem have been given a precise 
mathematical sense, and are computed by means of a uniquely defined prescription. Yet, the 
appearance of the undetermined constant A suggests that the present formalism might be 
physically incomplete, at least in this present stage of development. Notice that the constant 
A enters only the term proportional to G A m\m\(rrii + m 2 ) in the 3PN energy, and that for 
general orbits, the energy contains also 164 other terms which are all uniquely determined. 
The details of these calculations will be published elsewhere. 



(5.13) in Fl ). We find that there is perfect agreement provided that Static 



11 \ _ 1987 
3 840 



20| . Let us stress that in the present formalism we do not meet any ambiguity in the 
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